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UNIT  ONE 


INTRODUCTION  TO  MAPPING,  CHARTING,  AND  GEODESY  (MC&G) 


SUGGESTED  LESSON  PLANS* 


Lecture 

Material  Covered 

Textbook 

1 

Introduction  and  Orientation 

1.1 

2 

Ellipsoid  and  Geoid 

1.2.1 

3,  4,  5 

Surveying  and  Geodetic  Positioning 

1.2.2 

6,  7 

Photogramme try 

1.2.3 

8,  9 

Reference  Systems  (Earth-Fixed  and 
Space-Fixed);  Time 

1.2. 4.1 

1.2. 4. 2 

10 

Astronomic  Positioning 

1.2.4. 3 

11,  12 

Geodetic  Datums 

1.2.5 

13,  14 

Map  Projections 

1.2.6 

15 

Inertial  Survey  Systems;  Radio 
Interferometry 

1.2.7 

16,  17 

Gravity  and  Geoid 

1.3.1, 

1.3.2 

18 

Gravimetry  and  Gradiometry 

1.3.3 

19,  20 

Satellite  Geodesy 

1.4 

*Note ;  These  lesson  plans  are  based  on  42-lectures/semester, 
with  Units  One  and  Two  covered  in  the  first  semester; 
Units  Three  and  Four,  in  the  second. 
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UNIT  TWO 


Lecture 

21 

22  ,23 

24  ,25 

26,  27 
28 
29 

30-33 

34 

35-37 

38,  39 

40,  41 
42 


WEAPON  SYSTEMS  AND  MC&G 

SUGGESTED  LESSON  PLANS 

Material  Covered  Textbook 

Introduction  and  Orientation  2.1 

World  Geodetic  Systems  2.2;  also 

1.2.5 
(Unit  One) 

Introduction  of  Inertial  Navigation, 

Mechanization  of  Inertial  Navigation 


Systems  2.3,  2.3.1 

> 

INS  Errors  2.3.2 

Effects  of  Gravity  on  INS  2.3.3 

Gravity  Models  2.3.4; 

also  1.3.2 
(Unit  One) 

Initialization  of  Inertial  Systems  2.3.5 

Launch  Point  Errors  2.3.6 

Techniques  for  Aided  Inertial  Naviga¬ 
tion  Examples  2.3.7 

Military  Mapping  and  Charting  2.4 

Digital  Data  Bases  2.5 

Mission  Planning  •  2.6 
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Lecture 

UNIT  THREE 

SOURCE  DATA  COLLECTION  AND  REMOTE  SENSING 

SUGGESTED  LESSON  PLANS 

Material  Covered 

Textbook 

43 

Introduction  to  Remote  Sensing 

3.1 

44, 

45 

Photographic  Systems 

3.2.1 

46, 

47 

Electro-optical  Imaging  Systems 

3.2.2 

48 

Digital  Image  Characteristics; 
Non-Silver  Materials 

3.2.3, 

3.2.4 

49 

Other  Photographic  Systems;  Modula¬ 
tion  Transfer  Function 

3.2.5, 

3.2.6 

50 

Basic  Imaging  Operations 

3.3.1, 

3.3.2 

51 

Scanning 

3.3.3 

52, 

53 

Image  Enhancement 

3.3.4 

54, 

55 

Atmospheric  Effects 

3.4 

56 

Properties  of  Infrared  Radiation 

3.5.1 

57 

Infrared  Sensors 

3.5.2 

58, 

59 

Active  and  Passive  IR  Imagery 

3.5.3, 

3.5.4 

60 

HR  Scanning 

3.5.5 

61 

Properties  of  MW  and  MMW  Radiation 

3.6.1 

62, 

63 

Passive  Sensors  and  Imagery 

3.6.2, 

3.6.3 

64 

Radar  Sensors 

3.6.4 

65 

Active  Imagery;  SAR 

3.6.5 

66 

Target  Positioning 

3.7 
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UNIT  FOUR 


SATELLITE  SYSTEMS 


SUGGESTED  LESSON  PLANS 


Lecture 

Material  Covered 

Textbook 

67 

Introduction  to  Satellite  Systems 

4.1 

68, 

69, 

70 

Satellite  Platform  Considerations 

4.2 

71 

Gravitational  Perturbation  Effects 

4.3 

72, 

73 

Attitude  Reference  and  Control 

4.4.1,  4.4.2 
4.4.3 

74 

Satellite  Orbit  Control 

4.4.4,  4.4.5 

75 

Ground  Tracking  Systems 

4.5.1 

76 

Satellite  Altimetry 

4.5.2 

77 

Satellite-to-Satellite  Tracking 

4.5.3 

78, 

79, 

80 

Navy  Navigation  Satellite  System 

4.6 

81, 

82 

Global  Positioning  System 

4.7 

83, 

84 

Remote  Sensing  and  LANDSAT 

4.8 
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UNIT  ONE 

REVIEW  EXERCISES 

Chapter  Two 
(Section  1.2.1) 

1.  At  what  latitude  is  the  difference  between  geodetic  and 
geocentric  latitude  a  maximum? 

A:  A  =  ( (|> - *  )  =  (f  +  ^  f2)  sin  24>  -  j  f2  sin  4  <J>  (1-1) 

^  =  2 ( f  +  ~  f2)  cos  24.  -  2f2  cos  4«(.  (1-2) 

The  extremum  occurs  when  =  0. 

d<j> 

2 ( f  +  ~  f2)  cos  2<|.  =  2f2  cos  44.  (1-3) 

(1  +  ^  f)  cos  24>  =  f  cos  44>  (1-4) 

Let  x  =  cos  24> 

Then  cos  44>  =  2  cos2(24>)  -  1  =  2x2  -  1  (1-5) 

This  leads  to  the  quadratic  equation: 

(1+2  f)x  =  f[2x2  -  1]  (1-6) 
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f 


(1-7) 


(1  +  ^  f)x  =  2fx2  - 
(2f)x2  -  (1  +  \  f)x  -  f  =  0 

(1  +  \  f)  ±  ((1  +  \  f)2  -  A(2f)(-f)]^ 
x  =  __ 

(1  +  5  f)  ±  [l  +  f  +  |^  f2]^ 


(1-8) 


(1-9) 


(1-10) 


For  f  =  1/298.26  (WGS  72) 
x  =  149.38  (extraneous  root)  or  -0 . 013388 

If  x  =  cos  2<p  =  -0.013388,  then 


2«J)  =  90.767116 

(1-11) 

(1)  =  45.38 

(1-12) 

At  which  value 

A  =  0.19  deg 

(1-13) 

Using  the  facts  that  A 

=  0  at  latitudes  of  0  deg  and  90 

deg,  that  there  is  only  one  extremum,  and  that  A  is  po¬ 
sitive  at  the  extremum,  the  student  may  easily  demon¬ 
strate  that  this  extremum  is  a  maximum. 

Of  course,  the  student  may  prefer  to  solve  Eq. 
(1-2)  by  an  appropriate  numerical  method,  in  which  case 
the  extremum  is  verified  to  be  a  maximum  by  evaluation  of 
Eq .  (1-1)  at  points  in  the  vicinity  of  the  extremum. 


(Section  1.2.2) 


2.  A  tri angulation  survey  starts  with  a  baseline  (running 
west  to  east)  of  2000  m  connecting  points  Pq  and  .  At 
point  Pq ,  the  angle  between  the  baseline  and  an  unknown 
point,  P2 ,  is  measured  as  44  deg.  At  point  P-^  ,  the  angle 
between  the  baseline  and  point  ?2  is  measured  to  be  106 
deg.  Using  Pq  as  origin,  what  are  the  coordinates  of 
point  P2? 

A:  Using  the  formulas  for  plane  triangles,  it  is  found  that 

?2  lies  2671.00  m  north  of  Pq ,  2765.90  m  east  of  Pq .  The 
line  PqP^  *s  3845.05  m  in  length;  2778.63  m. 

3.  Continuing  the  triangulation  survey  of  Exercise  2,  the 
line  Pj?2  is  used  as  a  baseline  to  survey  a  new  point, 

P^.  At  P^ ,  the  angle  between  the  baseline  (P^^^nd 
point  P^  is  52  deg.  At  P2 ,  the  angle  between  the  base¬ 
line  and  P^  is  88  deg.  Find  the  coordinates  of  P^  rela¬ 
tive  to  the  origin,  Pq . 

A:  P^  is  6005.57  m  north  and  1618.36  m  east  of  Pq .  Line 

P2P3  is  3406.40  m;  line  ? 2  ^320.15  m  in  length. 

4.  A  trilateration  survey  starts  with  a  baseline  from  the 
origin,  Pq ,  to  an  initial  point,  P^ ,  with  length  900  m 
and  azimuth  (measured  positive  from  north  toward  east)  of 
30  deg.  The  following  distances  are  measured: 


FROM 

TO 

MEASURED  DISTANCE 
(m) 

Pi 

P2 

1000 

P0 

P2 

1100 

P2 

P3 

1200 

P1 

P3 

1800 

7 


Find  the  coordinates  of  »  an<*  ^3  relative  to  the 

origin  (PQ). 

A:  Using  standard  formulas  of  plane  trigonometry: 


DISTANCE  NORTH 

DISTANCE  EAST 

POINT 

op  po 

(m) 

op  po 

(m) 

P1 

779.42 

450.00 

P2 

19.39 

1099.83 

P3 

450.60 

2219.68 

5.  An  open  traverse  (refer  to  Fig.  1.2-9)  begins  at  the  ori 
gin,  PQ,  along  an  initial  baseline  with  azimuth  75  deg 
and  length  1600  m,  connecting  Pq  with  P^.  The  traverse 
continues  as  shown: 


FROM 

TO 

MEASURED  DISTANCE 
(m) 

MEASURED  ANGLE 

(deft) 

P1 

P2 

2000 

160 

P2 

P3 

800 

200 

P3 

P4 

2500 

250 

Angles  are  measured  in  a  clockwise  sense  from  the  previ¬ 
ous  traverse  line.  Find  the  coordinates  of  points  P-^ , 
?2>  P3>  and  P4  relative  to  the  origin  (Pq). 

A:  From  the  formulas  for  a  right  triangle: 
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POINT 

DISTANCE  NORTH 

op  po 

(m) 

DISTANCE  EAST 

OF  PQ 
(m)U 

P1 

414.11 

1545.48 

P2 

1561.26 

3183.79 

P3 

1768.32 

3956.53 

P4 

334.38 

6004.41 

6.  Referring  to  Exercise  5,  suppose  that  distance  measure¬ 
ments  are  correct  to  within  0.1  percent,  while  angles 
are  correct  to  within  1.0  min.  Find  upper  bounds  for  the 
position  error  of  point  . 


A:  The  upper  bound  error  analysis  is  based  on  the  formulas 

for  a  right  triangle: 


x  = 

s 

cos  e 

(6-1) 

and 

y  = 

s 

sin  9 

(6-2) 

To  first 

order , 

6x  = 

(cos  6 )6s 

- 

s  sin  6  60 

(6-3) 

l6x  |  < 

jcos  9  j !  6  s i 

+ 

[s  sin  e||6e| 

(6-4) 

6y  = 

(sin  9 )5s 

+ 

s  cos  0  60 

(6-5) 

ky  1  = 

|  sin  0 j [ 6 s  J 

+ 

js  cos  9 ) |<5 0 1 

(6-6) 

where : 


s  =  measured  distance 
9  =  measured  angle 

6s, 69  =  errors  in  measured  quantities 
6x,6y  =  errors  in  computed  coordinates 
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Applying  these  formulas  in  succession  to  the  four  right 
triangles  associated  with  the  traverse  gives  these  upper 
bounds : 


FROM 

TO 

NORTH  ERROR 
(m) 

EAST  ERROR 
(m) 

po 

P1 

0.86 

1.67 

P1 

P2 

1.62 

1.97 

P2 

P3 

0.43 

0.83 

P3 

P4 

2.03 

2.47 

P0 

P4 

4.94 

6.94 

For  a  statistical* approach  to  traverse  errors,  use  the 
methods  in  Chapter  Two  of  Introduction  to  Surveying  by 
Mueller  and  Ramsayer. 

(Section  1.2.3) 


7.  Describe  the  factors  that  limit  the  accuracy  of  maps  pre¬ 
pared  by  photogrammetric  techniques  from  overlapping  aerial 
photographs . 


A:  To  answer  this  question  completely,  the  student  should 

review  all  the  steps  in  the  process  of  map  preparation 
(Section  1.2.3).  Key  sources  of  error  that  the  student 
should  identify  are  those  associated  with 


•  The  camera  (lens  distortion,  film  plane 
flatness  and  stability,  etc.) 

•  Aircraft  position  and  attitude  control 

•  Horizontal  and  vertical  control  survey 
errors 

•  Operator  and  equipment  limitations  in  the 
stereocompilation  process. 
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(Section  1.2. A) 


8.  It  is  stated  in  the  text  (Section  1.2. 4.1)  that  "a  per¬ 
fect  atomic  clock  would  not,  after  a  lapse  of  many  years, 
correctly  predict  such  phenomena  as  sunrise,  sunset,  star 
transits,  eclipses,  etc."  Explain  why  this  is  the  case. 

A:  Because  these  phenomena  depend  on  earth  rotation,  which 

is  nonuniform. 

9.  Using  the  Star  Catalog  section  reproduced  as  Table  1.2-7 
(Section  1.2. 4. 2),  calculate  the  angular  distance  between 
the  stars  a  CMa  (Sirius)  and  y  Gem. 

A:  The  spherical  triangle  to  be  solved  has  sides  of 

90°  -  61  =  73.58 

and 

90°  -  «2  =  106.69 
with  an  included  angle  of 

<*2  -  =  1 .93 

The  third  side  is  found  most  easily  from  the  Law  of  Co¬ 
sines,  and  is 

d  =  33.16 
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If  an  individual  weighs  exactly  80  kg  at  sea  level  on 
the  equator,  how  much  would  this  person  weight  at  the 
North  Pole? 

A  correct  answer  requires  an  understanding  of  the  distinc¬ 
tion  between  weight  and  mass .  The  kilogram  is  a  unit  of 
mass,  and  an  80-kg  person  has  the  same  mass  anywhere.  The 
term  weight  is  a  measure  of  the  force  exerted  by  gravity: 

f  =  ma 

where : 

f  =  force  in  newtons  (N) 
m  =  mass  in  kg  2 

a  =  acceleration  in  gravity  in  m/sec 

2 

Using  9.7803  m/sec  for  the  acceleration  of  gravity  at 
the  equator  (Eq.  1.3-10,  Chapter  Three) 

f  =  782.42  N 

Since  scales  are  calibrated  in  equivalent  mass  units, 
rather  than  force  units,  the  force  of  782.42  N  is  indi¬ 
cated  as  80  kg.  At  the  pole  (Eq.  1.3-10),  the  accelera- 

2 

tion  of  gravity  is  9.8321  m/sec  ;  hence 
f  =  786.57  N 

and  the  scale  would  read 

m  =  80.42  kg 

The  apparent  excess  weight  (in  mass  units)  is  420  g 
(nearly  1  lb ) . 


11.  If  the  earth  were  to  stop  rotating,  by  how  much  would 
gravity  at  the  equator  change? 


A:  The  centrifugal  acceleration  is 


_  2 
a  -  iu 


r  cos 


with 

u>  =  7 . 3  x  10  5  rad/sec  (earth's  rotation  rate) 
r  =  6.378  x  10~*  m  (earth's  equatorial  radius) 

'  =  0  ( latitude ) 


Hence 

a  =  0.034  m/sec^ 

This  is  about  0.35  percent  of  the  nominal  sea  level  value. 

12.  Compute  the  excess  gravity  acceleration  (gravity  anomaly) 

at  the  surface  of  the  earth,  caused  by  a  sphere  of  depleted 
uranium,  with  a  mass  of  100  kg,  buried  10  m  below  the 
surface.  Would  the  presence  of  this  object  be  detectable 
by  the  use  of  ordinary  gravimeters? 

A:  Because  gravity  accelerations  are  additive,  the  excess 

acceleration  (anomaly)  directly  above  the  object  can  be 
computed  as  the  acceleration  due  to  the  buried  sphere 
alone,  ignoring  the  attraction  of  the  earth.  From  Ap¬ 
pendix  A  (Eq.  A. 1-6),  the  acceleration  is 


where : 
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G  =  6.67  x  10  ^  m'Vkg'sec^ 
M  =  100  kg 
h  =  10  m 


since  the  acceleration  due  to  a  uniform  sphere  is  the 
same  as  that  due  to  a  point  mass  located  at  the  center  of 
the  sphere.  (The  fact  that  the  sphere  is  uranium  is,  of 
course,  irrelevant.)  Therefore, 


..  _  (6.67xl0-11)(100)  _  .  ,,  ,„-ll  ,  2 

a  -  -  - - - -  -  *6.67  x  10  m/s ec 

(IQ)2 


Since  10  m/sec^  =  1  mgal, 


a  =  -6.67  x  10”  mgal 


The  excess  acceleration  would  not  be  measurable. 


13.  Referring  to  Appendix  B  as  required,  verify  Fig.  1.3-4b, 


A:  Using  Eq .  1.2-16  with  n=2;  m=2  gives 


^221  =  f>22^s:’'n  ^ ^  s:*"n  ^ 

=  3  cos^ '  sin  2A 


The  variation  along  a  meridian  (for  example,  A  =  45  deg) 
is  then  as  shown  in  Fig.  1.3-4b. 

14.  Referring  to  Appendix  B  as  required,  plot  the  Legendre 
polynomial  (sin  $')  and  verify  Fig.  1.3*4a. 
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A:  P^sin  0')  =  ^  sin^  0 '  -  ^  sin  0' 

=  j  sin  0'(5  sin^  0’  -  3) 

The  student  should  sketch  a  plot  of  this  function  between 
0’  =  -90  deg  and  <p '  =  90  deg,  showing  three  zero-crossings 
and  four  extrema.  This  plot  is  then  associated  with  the 
variation  of  the  spherical  harmonic  (n=3;  m=0)  along  a 
meridian . 

(Section  1.3.2) 

15.  What  is  the  normal  gravity  at  latitude  48  deg  N  using 
the  International  Gravity  Formula  [Eq.  (1.3-5)]? 

A.  980899.78  mgal 

16.  What  is  the  normal  gravity  at  latitude  48  deg  N  using 
the  Geodetic  Reference  System  1967  gravity  formula? 

A:  980890.15  mgal 

17.  At  a  point  at  sea  level  at  latitude  48  deg  N,  there 
is  a  measured  gravity  of  980933.8  mgal.  What  is  the 
gravity  anomaly  at  that  point? 

A:  Ag  =  g  -  y 

34.0  mgal  International  Gravity  Formula 
43.6  mgal  GRS  67 


( Section  1.3.3) 

18.  Suppose  that  a  gravity  measurement  is  to  be  made  using  a 
pendulum.  If  the  pendulum's  period  is  to  be  measured 
very  precisely,  how  accurately  must  the  dimensional  sta¬ 
bility  of  the  pendulum's  length  be  maintained  for  meas¬ 
urement  accuracy  of  one  mgal?  one  pgal? 


A:  Pendulum's  period  T  =  2njS./g.  Hence  g(measured)  = 


2 

4n  £ 


The  nominal  value  of  g  is  980  gal.  Thus  one  mgal  accu- 

-3 

racy  requires  10  /9S0^  one  part  in  one  million  dimen¬ 

sional  stability.  One  microgal  accuracy  requires  one 
part  in  a  billion. 

19.  In  Exercise  18,  what  is  the  timing  accuracy  required  to 
be  commensurate  with  a  one  mgal  level  of  dimensional 
stability . 


A:  Perturbing  the  equation 


g(measured)  = 


4n2£ 


gives 


.  /  JX  _  4n2A£  8n2JL  ^ 

Ag(measured)  =  — ~ —  -  — 5—  AT 


For  the  two  error  terms  to  be  of  the  same  magnitude,  it 
is  necessary  that 


4tt26£  _  8n2£AT 
t2  ~  t3 


Simplifying  gives 


A£  _  7  AT 
£  *  T 
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One  mgal  accuracy  requires  ££/£  to  be  one  part  per  mil¬ 
lion.  The  timing  uncertainty  for  one  mgal  of  timing- 
induced  error  is  therefore  one-half  part  per  million. 

20.  Suppose  that  a  pendulum  gravimeter  has  both  a  one  part 
per  million  (RMS)  dimensional  error  and  an  RMS  timing 
error  as  computed  in  the  previous  exercise.  If  both  of 
these  errors  are  random  and  are  unrelated,  what  is  the 
RMS  error  in  measured  gravity? 

A:  Independent  error  sources  add  by  summing  the  squares  of 

the  RMS  values  of  each  error  source.  The  square  root  of 
the  sum  is  the  RMS  output  error.  Both  length  and  timing 
errors  contribute  one  mgal.  Thus 

(1.02  +  l.O2)3*  =  1.4  mgal  RMS 
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UNIT  TWO 

REVIEW  EXERCISES 


A:  No,  any  convenient  reference  frame  is  suitable,  but 

if  it  is  non-inertial ,  the  body  forces  acting  on  the  ve¬ 
hicle  must  be  accounted  for  in  the  course  of  the  mechani¬ 
zation. 


2.  Why  is  it  not  possible  to  use  an  inertial  system's  accel¬ 
erometers  to  measure  and  subtract  away  the  effect  of 
gravity? 

A:  The  accelerometer  measurement  cannot  distinguish  be¬ 

tween  gravity  and  vehicle  acceleration.  This  is  a  funda¬ 
mental  physical  law  and  cannot  be  altered  by  improved 
accelerometers . 


3.  A  local ly- level  mechanized  inertial  system  provides  navi¬ 
gation  information  for  an  aircraft  flying  at  constant 
altitude  along  a  meridian  at  800  km/hr.  Neglecting  earth 
rotation,  what  is  the  angular  slew  rate  of  the  inertial 
platform  with  respect  to  inertial  space?  (Hint:  Con¬ 
sider  the  angle  between  the  plat  form- indicated  vertical 
and  the  earth's  spin  axis.) 
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A:  The  platform  maintains  itself  with  respect  to  the 

local  vertical  as  the  aircraft  moves.  Thus  the  platform's 
angular  rate  with  respect  to  the  earth's  center  is  the 
same  as  the  aircraft's  rate.  Since  the  motion  is  along  a 
meridian,  the  angle  between  the  earth's  spin  axis  and  the 
platform-indicated  vertical  changes  as 


ai rcraf  t 


R 


earth 


rad/hr 


3 

Using  6.4*10  km  as  an  approximate  value  for  the  equator¬ 
ial  radius  of  the  earth,  the  slew  rate  is 

•  =  800_km/hr_  =  013  rad/hr 
6. 4xl0'5  km 


4.  Consider  a  vehicle  with  a  local-level-mechanized  inertial 

system  which  is  stationary  at  the  equator,  at  a  place 

where  each  component  of  the  deflection  of  the  vertical  is 

10  sec  and  the  vertical  gravity  disturbance  is  50  mgal. 

Suppose  that  the  position  error  is  one  km  in  the  two  level 

channels  and  3.05  m  in  the  vertical.  Also  suppose  that 

the  gyro  errors  have  results  in  ij)  angle  errors  of  5  sec 

in  all  three  channels.  The  accelerometer  errors  are  5*10  ^ 

g  for  each  accelerometer  (where  the  g  is  a  unit  of  accel- 

2 

eration  equal  to  nominal  sea-level  gravity  --  9.8  m/sec  ). 
Tabulate  the  value  of  each  term  in  Eqs.  2.3-11  through 
2.3-12.  Note  that  these  terms  correspond  to  the  signals 
entering  the  acceleration  summing  nodes  in  Fig.  2.3-9. 

Use  units  of  10~  g.  What  does  this  example  illustrate 
about  the  cross-coupling  terms  of  Fig.  2.3-9? 
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Individual  Terms 

Accelerometer  Errors 
<mn.  m£.  MZ> 

5 

Cyro  Errors* 

<-8*e.  -g+N>  1 

i 

!  -95  (4.9  psec/sec! 

|  «PE.-  g*  Ch) 


Gravity  Field  Errors 
<°gN.  «*z) 


Velocity  Error 
(-20  sin  t.«VE, 

20  sin  l.6VN  +  20  cos  LSVj. 
-2  0  cos  L  6V„) 


VE  *  VZ  "  0 

2 

Suppose  an  inertial  system  uses  a  GM/r  gravity  model 
instead  of  a  reference  ellipsoid  gravity  model.  If  the 
nominal  surface  value  of  gravity  is  978.049  gal,  use  the 
approximate  ellipsoidal  gravity  formula. 

g  =  978.049(1  +  0.0053  sin2  <J> )  gal 

to  compute  the  vertical  acceleration  error  sensed  by  the 
system  at  a  latitude  of  45  deg.  At  this  rate  how  long 
does  it  take  for  100  m  of  altitude  error  to  accumulate? 
Assume  that  the  vertical  gravity  disturbance  is  zero  at 
the  location  of  the  system. 


A:  The  mechanized  value  of  gravity  is  978.049  cm/sec  . 

2 

The  ellipsoidal  value  of  g  at  45  deg  is  980.641  cm/sec  . 

2 

The  difference  between  the  two  is  2.6  cm/sec  .  Height 
12 

error  grows  as  ^  ( 2 . 6 ) t  .  One  hundred  meters  of  error  is 
reached  after  88  seconds. 


m 


6.  Suppose  a  guidance  system  is  to  be  compensated  for  gravi¬ 
ty  by  using  a  model  of  the  form  of  Eq.  2.3-20.  What  error 
sources  might  be  associated  with  such  a  model  if  the  maxi¬ 
mum  value  of  n  is  20. 

A:  Accuracy  of  the  value  of  GM ,  accuracy  of  C  and  S 

J  nm  nm 

coefficients  (errors  of  commission),  and  the  contribu¬ 
tions  of  the  neglected  terms  beyond  degree  20  (errors  of 
omission ) . 

(Section  2.3.5) 

7.  In  a  ballistic  missile  application  it  is  desired  to  erect 
the  inertial  guidance  platform  with  the  accelerometer 
sensing  axes  symmetrically  distributed  about  the  local 
vertical  (gravity  vector)  and  with  one  accelerometer  sens¬ 
ing  axis  pointing  downrange .  Assuming  the  use  of  a  mutual¬ 
ly  orthogonal  set  of  accelerometers,  what  are  the  nominal 

(error-free)  accelerometer  outputs  (f  ,  f  ,  and  f  )  when 

x  y  z 

alignment  is  complete?  (Hint:  Use  Eq .  2.3-23  and  symmetry 
considerations.  ) 

A:  f  =  f  =  f  =  0.58  g. 

x  y  z  & 

8.  A  local  level  inertial  platform  is  aligned  in  a  station¬ 
ary  vehicle  with  the  (orthogonal)  x  and  y  sensor  axes  in 
the  level  plane  at  a  point  on  the  earth's  surface  at  which 
the  geodetic  and  astronomic  verticals  coincide  (zero  de¬ 
flections  of  the  vertical)  and  the  common  latitude  is  <J> 
deg.  The  x  gyro  axis  points  5  deg  north  of  east.  The  gyro 
torquing  rate  commands  from  the  navigation  computer  are 
given  by: 


21 


0 


uj  =  ft  cos  0 

y 

uu  =  ft  sin  <Ji 
z 

where  ft  =  Earth  Rotation  Rate  relative  to  Inertial  Space. 

What  are  the  observed  platform  tilt  rates  (<i>  ,  and  <P  ) 

*  > 

about  the  platform  x  and  y  axes  respectively?  (Hint: 
Compute  the  components  of  earth  rate  appearing  about  the 
actual  (misaligned)  sensor  axes  (ft^  and  fty^)  and  use  Eq  • 
2.3-27  to  get  tilt  rates.) 

A:  <i>  =  ft  cos  <p  sin  5° 

x 

<i>  =  ft  cos  <j>  ( 1  -  cos  5°) 

y 

Using  the  small-angle  approxima t ions ,  sin  5°  a  0.1  and 
cos  5°  «  1.0,  evaluare  the  tilt  rates  <i>x  and  <i>^  derived 
in  Problem  8  for  vehicle  latitudes  =  30  deg  N  and  (j>2  = 
70  deg  N.  Answer  the  following  questions: 

(1)  Which  accelerometer  senses  a  tilt  rate 
that  can  be  used  to  estimate  and/or  con¬ 
trol  the  platform  azimuth  misalignment? 

(2)  What  is  the  change  in  the  magnitude  of 
the  gy rocompassing  control  signal  as  the 
alignment  latitude  increases? 

( 3 >  What  does  this  difference  tell  you  about 
inertial  platform  sel f-al ignment  capabil¬ 
ities  at  high  latitudes  in  the  presence 
of  accelerometer  time-varying  errors  or 
gyro  bias  drift  errors? 
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A: 

At 

latitude  <)>  ^  : 

4>  = 

X 

0.087  0  and  4>  =0 

y 

At 

latitude  q>2: 

i  = 

X 

0 . 034  fi  and  <i>  =0 

y 

(1) 

The 

y- Accelerometer . 

Tilt  rates  appear  about 

the 

x-axis  only. 

(2) 

The 

control  signal  (4> 

)  decreases  with  increas- 

ing 

latitude . 

(3) 

The 

decrease  in 

the  control  signal  with  increas 

ing 

latitude  implies 

increasing  difficulty  in 

reaching  a  satisfactory  azimuth  alignment 
solution  in  the  face  of  accelerometer  and 
gyro  errors.  The  "signal- to-noise”  ratio 
decreases  at  higher  latitudes.  This  is  a 
manifestation  of  approaching  collinearity 
between  earth  spin  rate  and  local  gravity 
vectors . 


10.  An  error-free,  north- slaved ,  local  level  inertial  platform 
is  self-aligned  in  a  stationary  vehicle  at  a  point  on  the 
earth's  surface  at  latitude  30  deg  N.  An  easterly  deflec¬ 
tion  of  the  vertical,  n  =10  sec,  exists  at  the  alignment 
location.  The  existence  of  this  deflection  of  the  vertical 
is  not  accounted  for  in  the  navigation  computer.  Evaluate 
the  effect  of  this  omission  on: 


(1)  Platform  tilt  errors 

(2)  Platform  azimuth  error 

relative  to  the  (geodetic)  navigation  coordinate  frame  when 
alignment  is  complete.  Repeat  the  evaluation  when  the  ve¬ 
hicle  latitude  is  70  deg  N  and  the  same  conditions  exist. 
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A:  (a)  <t>  A  =  10  sec;  $  A  =  0 ;  Azimuth  error,  ii>  A  =  5.8  sec 

yO  xO  zO 

(b)  <S>yQ  =  10  sec;  =  0;  Azimuth  error,  =  27.5  sec 

Equations  2.3-28  gives  azimuth  errors. 

An  aircraft  is  equipped  with  an  avionics  suite  that  in¬ 
cludes  a  gimballed  inertial  platform,  a  Doppler  radar  set 
and  a  LORAN  receiver.  Two  distinct  and  separate  modes  of 
in-flight  initialization  are  mechanized  in  the  system: 

(1)  Doppler  (Velocity)  -  Inertial 

(2)  LORAN  (Position)  -  Inertial. 


In  the  former  mechanization  the  Doppler  velocity  inputs 
are  transformed  directly  into  inertial  sensor  coordinates, 
using  resolvers  mounted  between  the  gimbals  of  the  iner¬ 
tial  platform,  prior  to  comparison  with  the  inertial  sys¬ 
tem  velocity  outputs.  In  the  latter  mechanization  the 
LORAN  receiver  outputs  are  converted  to  geodetic  latitude 
and  longitude  in  the  navigation  computer  and  compared 
with  the  latitude  and  longitude  outputs  of  the  inertial 
system  to  form  alignment  control  signals. 

(a)  What  is  the  basic  difference  in  the  dy¬ 
namics  of  the  inertial  platform  alignment 
process  between  the  two  mechanizations? 

(b)  Can  a  complete  inertial  system  initializa¬ 
tion  be  achieved  if  the  LORAN  receiver 

is  inoperative? 

A:  (a)  In  the  second  mechanization,  azimuth  alignment 

varies  with  aircraft  flight  path  (spatial  rate 
gyrocompass i ng ) . 
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(b)  Not  unless  an  additional  method  for  discrete  or 
continuous  position  updating  is  available. 

(Section  2.3.6) 

12.  It  is  1985!  USAF  is  considering  the  deployment  of  a  new 
ballistic  missile  system  as  an  augmentation  to  the  aging 
Minuteman  force.  In  order  to  obstruct  a  surprise  attack 
by  opposing  forces,  this  new  missile  will  be  deployed  in 
truck-like  transports  which  will  roam  at  random  through 
the  deserts  of  southwestern  CONUS.  The  missiles  will  be 
erected  and  fired  from  the  moving  transporter  when  the 
order  to  launch  is  received. 

A  weapon  system  CEP  of  75  m  is  desired,  and  an 
allowance  for  the  guidance  system  contribution  to  this 
total  has  been  set  at  30  m,  including  the  effects  of  launch 
site  uncertainties.  The  latest  generation  of  inertial 
guidance  hardware  indicates  that  there  will  be  no  problem 
in  meeting  this  accuracy  goal  if  the  initial  condition 
errors  can  be  controlled.  The  transporter  has  according¬ 
ly  been  equipped  with  special  radio  devices  that  allow  it 
to  determine  its  position  continuously  to  an  accuracy  of 
about  3  m.  Has  the  initial  condition  problem  been  fully 
addressed?  Discuss. 

A:  No.  The  miss  sensitivities  listed  in  Table  2.3-1 

show  that  suppressing  the  launch  point  position  error  to 

a  level  of  3  m  is  not  sufficient.  The  velocity  error 

_  2 

also  will  have  to  be  constrained  to  the  order  of  10 
m/sec  if  the  accuracy  goal  is  to  be  achieved! 
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( Sec  t ion  2.3.7) 


13.  Why  are  aided  inertial  navigation  techniques  preferred 
over  unaided  inertial  navigation  techniques?  Give  exam¬ 
ples  of  the  two  types  of  Aided  INS  errors,  and  how  they 
are  counteracted. 

A:  Errors  in  an  unaided  INS  typically  grow  with  time  in 

an  unbounded  manner.  To  maintain  accuracy,  the  INS  used 
in  modern  aircraft,  missiles,  and  submarines  are  normally 
aided  or  "damped"  with  data  from  external  (i.e.,  non- 
inertial)  sources. 

The  two  types  of  Aided  INS  errors  are  those  1) 
associated  with  the  measurement  device  itself  and  referred 
to  as  instrument  errors  and  those  2)  induced  by  external 
factors  which  operate  independently  of  instrument  quality 
and  referred  to  as  environmental  errors. 

Examples  of  1)  are  gyro  drift,  accelerometer 
bias,  or  radio  receiver  noise. 

Example  of  2)  are  radio-wave  propagation  anoma¬ 
lies  and  barometric  variations  due  to  weather  patterns. 

14.  What  is  the  fundamental  technique  behind  terrain-matching 
navigational  methods?  Give  examples. 

A:  Terrain  matching  is  a  navigational  process  based  on 

recognition  of  terrain  features.  This  is  fundamentally  a 
differencing  technique,  where  position  information  is 
inferred  by  comparing  certain  measured  terrain  features 
(or  possibly  man-made  objects)  with  reference  features 
stored  in  some  type  of  map. 


Zo 


One  example  is  a  pilot  flying  under  VFR  who  es¬ 
tablishes  his  position  by  comparing  observed  ground  fea¬ 
tures  with  those  indicated  on  his  navigator's  map. 

Another  example  is  photographic  correlation 
schemes  of  map  matching.  With  these  schemes,  the  terrain 
is  observed  photographically  and  compared  with  photographi¬ 
cally-based  prestored  reference  maps  until  some  best  cor¬ 
relation  of  images  is  achieved  between  the  observation 
and  the  reference. 

15.  Autonomous  terminal-homing  guidance  systems  generally  use 
scene-matching  concepts  to  determine  the  location  of  the 
vehicle  or  its  position  relative  to  the  target.  Give 
three  specific,  different  difficulties  that  can  be  en¬ 
countered  in  preparing  and  using  the  reference  imagery 
model  by  the  vehicle’s  onboard  navigation  system  for  cor¬ 
relation  guidance. 

A:  The  most  common  form  of  course  data,  from  which  each 

mission-dependent  reference  map  is  prepared,  is  stereo¬ 
scopic  optical  photography  from  high-altitude  aircraft  or 
satellites . 

One  problem  is  that  the  original  scene  viewed  in 
the  initial  photography  may  have  to  be  regenerated  as  if 
viewed  by  a  different  sensor  operating  at  a  different 
wavelength,  in  order  to  match  the  image  expected  from  the 
sensor  aboard  the  autonomously-guided  vehicle. 

Another  problem  is  the  difficulty  in  predicting 
imagery  at  the  same  wavelength.  For  example,  imagery  of 
a  building  complex,  obtained  with  a  high- resolut ion ,  ther¬ 
mal  IR  scanner,  is  dependent  in  part  on  such  unknown  factors 
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as  interna]  building  structure,  heating/air  conditioning 
of  the  buildings  during  the  period  prior  to  imaging,  re¬ 
cent  solar  illumination  history,  cloud  cover,  surface 
conditions,  fog/haze/smoke  densities,  etc. 

A  third  problem  is  that  the  scene  correlation 
algorithms  must  be  insensitive  to  image  distortion  aris¬ 
ing  from  variations  in  perspective. 

16.  Suppose  that  a  local-level  mechanized  inertial  platform 
points  toward  a  star  located  at  the  zenith.  Discuss  the 
improvement  in  system  azimuth  error  that  can  be  gained 
from  this  measurement. 

A:  Refer  to  Fig.  2.3-24.  The  star  tracker's  barrel  axis 

is  vertical,  thus  the  measurement  provides  no  angle 
azimuth  information. 

17.  Consider  an  inertial  system  which  has  no  platform  mis¬ 
alignment  (i.e.,  6  =  0  in  Fig.  2.3-24).  Can  any  error- 
free  star  fix  correct  all  other  system  errors?  Discuss. 

A:  The  star  fix  provides  angle  information  in  only  two 

coordinates.  The  component  of  along  the  line  of  sight 
to  the  star  is  not  improved. 

Chapter  Four 

18.  How  can  Point  Positioning  Data  Bases  (PPDBs)  be  used  to 
support  tactical  bombing  operations? 

A:  Targets  identified  in  reconnaissance  photographs  can 

often  be  located  approximately  with  respect  to  nearby 
landmarks.  The  positions  of  the  landmarks  in  an  appro¬ 
priate  coordinate  system  can  be  determined  using  the  PPDB 


photochips  and  the  APPS .  The  target  location  can  then  be 
determined  in  the  chosen  coordinate  system. 


Chapter  5 

19.  What  are  the  key  considerations  for  selection  of  a  navi¬ 
gation/guidance  concept  based  on  an  on-board  reference  map? 

A:  (a)  Available  on-board  computing  power  and  memory 

volume 

(b)  Sensor  type(s)  and  resolution 

(c)  Reference  data  representation  selection. 

20.  What  are  the  principal  sources  of  navigation  errors  which 
arise  from  the  use  of  an  on-board  reference  map? 

A:  (a)  Map  source  material  resolution,  collection 

geometry,  and  ground  control 


(b)  Translation  from  source  to  reference  map  data,  es¬ 
pecially  if  source  is  based  on  a  different  sensor 


(c)  Temporal  variations  in  sensed  data  for  the  given 
terrain . 

21.  (a)  What  are  the  principal  tasks  to  be  accomplished  in 

preparing  reference  maps  for  specific  missions? 

A:  (a)  Scene  content  identification 

(b)  Reference  effectiveness  prediction 

(c)  Scene  selection. 
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(b)  What,  implications  do  these  have  for  the  tactical  or 
strategic  mission  planner? 


A:  ( a ) 

(b) 


Process  is  labor-intensive/special  skills  required 
Extensive  computer  support  required 


(c)  Over-specification  may  result  in  slow  response, 
overload  of  on-board  system 

(d)  Error  budget  must  be  kept  in  mind. 


22.  What  are  the  elements  of  an  observable  feature  in  a  ter¬ 
rain  data  base? 


A:  (a)  Lack  of  ambiguity  relative  to  neighboring  features 


(b)  Minimal  measurement  uncertainty  with  respect  to 
on-board  sensor. 

Chapter  Six 


23.  Give  the  three  important  interdependent  performance  cri¬ 
teria,  as  related  to  route  analysis  for  low-altitude, 
high-speed  weapon  delivery  vehicles  such  as  cruise  mis¬ 
siles  penetrating  enemy-controlled  areas. 

A:  The  important  performance  criteria  are  1)  probability 

of  detection  (being  neutralized  by  the  enemy),  2)  range 
loss,  and  3)  probability  of  crashing  (clobber). 


24. 


Refer  to  Eqs. 
cruise  missile 
in  achieved  cl 


2.6-5  and  2.6-6.  Assume  that  a  particular 
flies  with  a  standard  deviation  of  error 
earance  altitude  of  25  ft,  and  a 


JO 


standard  deviation  of  the  clearance  rate 

sec.  What  is  the  lowest  mean  flight  altitude  that  can  be 
achieved,  for  a  2-hour  flight,  that  results  in  at  most  a 
1%  probability  of  clobber? 

A: 


P  ,  ,n  =  0.01  =  1  -  e'AT  =  1  -  e'7200X 

clobbers  1 ) 

1 ,39588xl0*6/sec  =  A 

1 . 39588*  10*  ^  =  A  =  i—  exp - 

2n  °ACL  2o2ACL 


10.035  = 


2 


2o 


ACL 


Mean  Flight  AltitudeA  =  112  Feet 

M 
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UNIT  THREE 
REVIEW  EXERCISES 


Chapter  One 


1.  Define  the  term  remote  sensing  and  give  several  examples. 

A:  As  defined  in  Section  3.1.1,  remote  sensing  refers  to 

techniques  for  obtaining  information  about  objects,  areas, 
or  phenomena  by  analyzing  data  collected  by  devices  lo¬ 
cated  at  a  distance  from  the  object  of  investigation.  In 
the  examples  discussed  in  the  text,  remote  sensing  is 
based  on  electromagnetic  energy  that  is  either  emitted  or 
reflected  by  the  objects  under  study.  This  energy  may  be 
visible  light,  ultraviolet,  infrared,  or  microwave. 

Chapter  Two 

2.  Briefly  summarize  the  principal  characteristics,  advan¬ 
tages,  and  disadvantages  of  the  two  major  imaging  tech¬ 
nologies  --  photochemical  and  photo-electronic. 

A:  The  student's  answer  should  include  the  information 

contained  in  Table  3.2-1,  as  well  as  in  the  introduction 
to  Chapter  Two. 

3.  It  is  frequently  stated  that  fluorescent  light  sources 
are  more  efficient  than  incandescent  devices.  Explain 
the  meaning  of  the  statement. 
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A:  Efficiency,  in  this  context,  is  measured  by  the  lumen- 

to-watt  ratio  of  the  device.  For  fluorescent  lights  the 
ratio  can  be  as  high  as  80  lumens/watt,  while  a  typical 
figure  for  incandescent  light  is  10  lumens/watt.  The 
incandescent  lamp  dissipates  about  20  percent  of  the  energy 
supplied  to  it  in  the  form  of  heat  conducted  through  the 
socket  or  convected  into  the  surrounding  air.  The  remain¬ 
ing  80  percent  is  radiated  in  the  form  of  electromagnetic 
energy,  but  most  of  this  energy  (about  85  percent)  is  in 
the  infrared.  The  fluorescent  bulb  and  its  associated 
circuitry  produce  considerably  less  heat;  also,  the  spec¬ 
tral  distribution  of  the  light  emitted  by  the  bulb  is 
better  matched  to  the  relative  spectral  efficiency  of  the 
human  eye  (Fig.  3.2-2).  Nonetheless,  compared  to  the 
theoretical  maximum  of  686  lumens/watt,  neither  type  of 
device  would  be  considered  efficient . 

U.  What  would  be  the  color  of  the  light  from  a  hypothetical 
device  that  produces  light  from  electricity  at  close  to 
the  theoretical  optimum  efficiency  of  686  lumens/watt? 

A:  Such  a  hypothetical  device  would  emit  monochromatic 

light  at  the  wavelength  of  maximum  spectral  efficiency 
for  the  human  eye,  about  550  to  560  nm  (green).  Refer  to 
Fig.  3.2-2. 

5.  On  a  very  bright  day  when  the  sky  is  clear,  the  illuminance 
due  to  sunlight  may  be  as  high  as  10^  footcandle.  Express 
this  value  in  metric  units. 

A:  The  correspondence  is  given  in  Table  3.2-3.  The  illumi- 

5  1 

nance  due  to  sunlight  would  be  about  10  lux  (lumen/m"). 
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f>.  Although  silver  fluoride  belongs  to  the  class  of  compounds 
called  silver  halides  ,  it  is  not  mentioned  in  most  discus¬ 
sions  of  silver-based  photographic  technology.  Explain. 

A:  Because  silver  fluoride  is  highly  soluble  in  water, 

it  is  not  usable  in  most  photographic  applications. 

7.  What  is  a  latent  image  in  photography? 

A:  The  student's  reply  should  summarize  the  relevant 

discussion  in  Section  3.2.1.  The  key  sentence  (page  3-19): 
"The  distribution,  within  the  emulsion,  of  altered  and 
unaltered  halide  microcrystals ,  in  accordance  with  the 
light  intensities  of  the  image,  is  called  the  latent  image." 

8.  Assume  that  it  is  necessary  to  examine  an  exposed,  but 
undeveloped,  strip  of  aerial  film  to  obtain  information 
about  the  latent  images  contained  thereon.  What  are  some 
ways  in  which  this  can  be  done? 

A:  It  is  expected  that  the  student  will  state  unequivocal¬ 

ly  that  a  latent  image  can  be  detected  only  by  developing  it. 

9.  For  the  H  and  D  curve  shown  in  Fig.  3.2-5,  estimate  the 
value  of  gamma . 

A:  As  the  curve  is  drawn,  gamma  is  about  1.33,  somewhat 

higher  than  the  usual  range  of  values  for  pictorial  or 
aerial  photography. 

10.  F.xplain  the  reciprocity  law  and  define  the  term  rec  iproc  i  ty 
failure. 
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A:  The  reciprocity  law  is  an  approxima t ion  which  states 

that  the  photographic  effect  (as  measured  by  the  dens i ty 
resulting  from  standard  development)  is  a  function  of  the 
product  of  intensi ty  and  exposure  time  ( lux- seconds ) ,  and 
does  not  depend  on  the  individual  values  of  these  two 
quantities,  so  long  as  the  product  is  the  same.  For  most 
pictorial  and  aerial  films,  there  is  a  fairly  wide  range 
of  illuminance  levels  for  which  the  approximation  gives 
satisfactory  accuracy.  At  very  low  or  very  high  illumi¬ 
nance  levels,  the  film  characteristics  are  not  described 
well  by  the  reciprocity  law.  This  phenomenon  is  called 
reciprocity  failure. 

11.  Why  is  reciprocity  failure  a  more  serious  problem  with 
color  film  than  with  black  and  white  materials? 

A:  The  three  emulsion  layers  in  color  film  typically 

have  different  sensitivity  and  reciprocity  characteristics. 
Consequently,  reciprocity  failure  affects  color  balance. 

The  resulting  color  shifts  can  sometimes  be  compensated 
for  by  the  use  of  filters  (if  the  illuminance  deviates 
only  slightly  from  the  design  range),  but  severe  reci¬ 
procity  failure  generally  leads  to  unusable  results. 

12.  Distinguish  between  photoconductive  and  photoemi ssi ve 
devices . 

A:  Both  types  of  devices  convert  light  signals  to  electri¬ 

cal  signals,  but  they  exploit  different  physical  principles. 
In  a  photoconduc t ive  device,  the  light  energy  decreases 
the  electrical  resistance  of  a  suitable  semiconductor 
material.  Thus  the  current  in  a  circuit  including  such  a 
device  is  modu 1  a t ed  in  accordance  with  changes  in  light 
intensity.  Fhutoemi  sstve  devices,  on  the  other  hand. 
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depend  on  the  photoelectric  effect,  in  which  incident 
light  causes  the  ejection  of  electrons  from  the  surface 
of  the  photosensitive  material.  These  photoelectrons  are 
collected  to  produce  an  electrical  signal. 

13.  As  a  general  rule,  there  is  one  characteristic  that  dis¬ 
tinguishes  all  of  the  non-silver  photographic  processes 
from  conventional  silver-based  technology.  Explain. 

A:  The  distinguishing  characteristic  is  photographic 

speed ,  or  sensitivity  to  light.  Conventional  materials 
can  be  made  far  more  sensitive  to  light  than  any  non¬ 
silver  material . 

14.  Explain  what  is  meant  by  the  modulation  transfer  function 
( MTF ) . 


A:  The  student  may  answer  in  informal  terms  (pages  3-73 

to  3-75)  involving  the  degradation  of  the  image  of  a  bar 
pattern  test  object;  or  in  formal  terms  (pages  3-76  to 
3-82),  as  the  absolute  value  of  the  normalized  Fourier 
transform  of  the  point  spread  function  (PSF). 


Chapter  Three 


Distinguish  between  the  concepts 


A:  The  process  of  sampling  breaks  up  a  continuous  picture 

into  a  finite  number  of  discrete  picture  elements  (pixels). 
Each  pixel  is  characterized  by  an  average  value  of  a  physi¬ 
cal  variable  like  light  intensity  or  optical  density.  The 
process  of  digi t izing  (or  quantizing)  divides  the  range 
of  values  of  the  physical  variable  into  a  discrete  number 
of  steps  (the  number  of  steps  is  usually  chosen  to  be  a 
power  of  two ) . 
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16.  The  image  enhancement  techniques  of  Chapter  Three  are 

described  in  terms  of  sampled  and  digitized  images  avail¬ 
able  in  digital  form  (for  computer  storage  and  manipula¬ 
tion'*.  These  techniques  are,  therefore,  directly  appli¬ 
cable  to  imagery  produced  by  elec tro-optical  devices. 

How  cs  they  be  applied  to  images  in  continuous  form  on 
photographic  film'1 

A:  By  sampling,  digitizing,  and  storing  the  film  image, 

using  the  methods  of  Sections  3.3.2  and  3.3.3.  It  may  be 
of  interest  that  there  are  certain  techniques  that  can  be 
applied  directly  to  film  images  without  the  use  of  a  com¬ 
puter.  These  methods  are  purely  photochemical  in  nature, 
and  are  most  effective  for  gray  scale  manipulation,  con¬ 
trast  enhancement,  and  edge  sharpening.  While  they  do 
not  require  the  use  of  scanning,  digitizing,  and  computing 
equipment,  they  do  require  highly  trained  personnel  and 
darkroom  facilities. 


Chapter  Four 


17.  This  exercise  refers  to  Figs.  3.4-la,  b  and  3.4-2. 

each  of  the  environmental/weather  conditions  listed 
and  for  a  fixed  aperture  size,  specify  which  of  the 
pheric  windows  allows  <_  2  dB/km  total  one-way 
while  maximizing  the  angular  resolution 


Atmospheric  Windows:  a . 

b. 

c . 

d. 
e  . 
f . 

£• 

h. 


For 
below , 
atmos- 
attenuation 
obtainable . 

Visible  light 
Mid-IR  (3-5  pm) 

Far- 1R  (8-14  pm) 

220  GHz 
140  GHz 
94  GHz 
35  GHz 
<  20  GHz 
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Weather  Conditions: 


1.  Fog,  400-ro  visibility 

2.  Haze,  1.2-km  visibility 

3.  Fog,  100-m  visibility 

4.  Moderate  rain  (0.4  cm/hr) 

5.  Haze,  23-km  visibility 

6.  Very  heavy  rain  (5  cm/hr) 

A:  The  appropriate  atmospheric  windows  are  as 

follows : 

Environmental  Condition  Window 


1 

2 

3 

4 

5 

6 


d 

c 

e  (marginally) 

g 

b 

h 


18.  What  is  meant  by  Rayleigh  scattering?  Give  an  example. 

A:  Rayleigh  scattering  refers  to  the  scattering  of  elec¬ 

tromagnetic  energy  by  particles  which  are  small  with  re¬ 
spect  to  the  wavelength  of  the  scattered  energy.  The 
amount  of  scattering  is  proportional  to  the  fourth  power 
of  the  frequency.  An  example  is  the  scattering  of  visible 
solar  radiation  by  air  molecules  (leading  to  the  reddish 
appearance  of  the  sun  at  rising  and  setting,  and  to  the 
blueness  of  the  sky). 
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Chapter  Five 


19.  This  exercise  refers  to  Figs.  3.5-1  and  3.5-3,  and  involves 
two  objects: 

•  Hot  metal  housing  surrounding  a  tank 
exhaust  vent;  temperature,  700  K;  IR 
reflectivity,  0.4 

•  Black  asphalt  road;  temperature,  270  K; 

IR  reflectivity,  0.05 

For  each  material,  find  the  wavelength  at  which  the  peak 
of  its  thermal  radiation  curve  occurs.  Find  the  total 
radiated  energy  (watt/m  )  at  the  wavelength  of  peak  emis¬ 
sion.  Select  a  detector  material  that  would  be  most  suit¬ 
able  for  sensing  the  radiation  from  each  of  the  objects. 


A:  The  results  are  summarized  as  follows: 


OBJECT 

PEAK 

WAVELENGTH 

PEAK  ENERGY 

DETECTOR 

MATERIAL 

Hot  metal 

4.3  pm 

2.6*10^  watt/m^ 

Ge:Au  or  HgCdTe 

Asphalt  road 

11.1  pm 

3.5*10^  watt/m^ 

HgCdTe 

Chapter  Six 

20.  This  exercise  refers  to  a  spaceborne  millimeter-wave  sen¬ 
sor  (140  GHz)  for  passive  detection  of  atmospheric  vehi¬ 
cles  from  radiometric  contrast  temperature  measurements. 

a.  Assuming  an  orbital  altitude  of  800  km, 

compute  the  diameter  of  the  antenna  aper¬ 
ture  required  for  resolution  of  5-meter 
target  objects. 
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b.  If  the  sensor  were  replaced  by  one  operat¬ 
ing  at  10  pm,  what  would  be  the  required 
antenna  aperture  for  resolution  of  a 
5-meter  object? 

A:  a.  The  angular  resolution,  in  radians,  is  given  by 

A/D.  At  the  stated  frequency  of  140  GHz,  it  would  be 
necessary  to  have 

D  =  341  m 

to  satisfy  the  stated  resolution  requirements.  In  effect, 
the  result  is  that  a  millimeter-wave  sensor  at  140  GHz 
could  not  resolve  objects  as  small  as  5  m. 

b.  If  A  is  changed  to  10  pm,  then  the  required  resolution 
is  obtained  for 

D  =  1 . 6  m 

a  requirement  that  could  be  met  in  practice. 


40 


UNIT  FOUR 


REVIEW  EXERCISES 


Chapter  Two 

1.  Explain  what  is  meant  by  the  term  telemetry  or  telemetry 
data . 

Answer :  Telemetry  refers  to  the  transmission  of  instrument 

readings  (measurements)  of  all  kinds  --  for  example,  tempera¬ 
ture,  pressure,  or  acceleration  --  to  an  observer  located  at  a 
remote  site. 


2.  Assume  that  a  spacecraft  has  the  form  of  a  circular  cylin¬ 
der  1  m  in  diameter  and  iOO  m  long.  The  spacecraft  is  aligned 
with  the  axis  of  the  cylinder  pointing  toward  the  center  of 
the  earth,  and  is  in  a  circular  orbit  at  a  height  of  600  km. 
Compute  the  difference  in  gravitational  attraction  between  the 
two  ends  of  the  cylinder. 


,->wer :  Assuming  an  inverse  square  gravi  tat  ional  law,  the 

difference  in  gravitational  attraction  is 


6g 


GM 

(R+h)* 


GM 

(R+h+tr 


where 


GM  is  the  geocentric  gravitational  constant 
(Table  1.3-2,  Unit  Two) 

R  is  the  radius  of  the  earth 
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h  is  the  satellite's  height  above  the  earth 
£  is  the  length  of  the  satellite 


For  simplicity,  substitute 

p  =  R+h 


and  use  the  binomial  theorem  to  obtain  a  series  expansion  for 
6g  in  terms  of  the  small  quantity  x.  The  result  is 

og  -  — y~  +  terms  of  higher  order 
P 

Using  the  approximate  numerical  values 

GM  =  4.0  x  1014  m3/sec2 

p  =  R+h  =  7  x  10^  m 

x  =  -  =  1.4  x  10"5 
P 


the  answer  is 

6g  =  2.3  x  10  4  m/sec2  =  23  mgal 

3.  If  an  earth  satellite  has  no  internal  heat  sources,  will 
it  eventually  cool  to  a  temperature  of  absolute  zero? 

Answer :  No.  It  will  reach  an  equilibrium  temperature  at  which 

the  solar  energy  absorbed  by  the  spacecraft  equals  the  infrared 
energy  emi tied  by  the  satellite.  The  value  of  the  equilibrium 
temperature  depends  on  the  reflective  and  absorptive  qualities 
of  the  satellite's  surface,  but  will  probably  be  in  the  range 
of  200  to  250  K  for  a  satellite  near  the  earth. 
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Chapter  Three 


4.  Assume  a  satellite  in  a  retrograde,  nearly  polar,  circular 
orbit  with  an  inclination  of  99  deg.  The  height  of  the  satellite 
at  its  equatorial  crossings  is  890  km.  Compute  the  daily  change 
in  Q,  the  right  ascension  of  the  ascending  node. 

An swer :  Using  6378160  m  for  the  earth's  equatorial  radius, 
the  satellite's  semi-major  axis  is  1.1395.  From  Eq .  (4.3-17) 
tne  node  advances  by  0.9864  deg/day. 

5.  For  the  satellite  described  in  the  preceding  exercise, 
what  is  the  annual  nodal  change?  Why  is  such  a  satellite  called 
sun -synchronous?  Why  are  sun- synchronous  satellites  important 
for  photographic  and  remote-sensing  applications? 

Answer:  The  annual  change  in  Q  is  very  nearly  360  deg.  The 

plane  of  the  satellite's  orbit  rotates  (with  respect  to  a  geo¬ 
centric  inertial  coordinate  system)  at  the  same  rate  as  the 
apparent  motion  of  the  sun  --  hence  the  name  sun-synchronous . 

In  practice,  this  means  that  if  the  satellite  crosses  the 
equator  on  a  northbound  pass  at  local  noon,  it  will  make  every 
northbound  pass  at  local  noon.  (Southbound  passes  will  be 
made,  of  course,  at  local  midnight.)  The  significance  for 
photographic  and  remote-sensing  applications  is  that  the  satel¬ 
lite  observes  the  earth,  at  every  pass,  under  virtually  iden¬ 
tical  conditions  of  solar  illumination. 

6.  Could  a  sun- synchronous  system  be  designed  involving  a 
satellite  in  a  direct ,  rather  than  retrograde ,  orbit? 

Answer:  No.  Examination  of  Eq .  (4.3-17)  shows  that  the  re¬ 

quired  posi ti ve  nodal  rate  can  be  obtained  only  if  cos  i  is 
negative,  which  implies  a  retrograde  orbit  ( i  >  90  deg). 
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Chapter  Four 


7.  What  are  some  of  the  means  of  attitude  determination  used 
by  satellites? 


Answer:  The  student's  answer  should  include: 


•  Sun  sensors 

•  Horizon  sensors 

•  Magnetometers  (to  detect  spacecraft  orien¬ 

tation  relative  to  the  earth's  magnetic 
field) 

•  Star  sensors 

•  Inertial  reference  systems. 


8.  Could  a  star  tracker  be  used  to  track  planets  instead  of 
stars?  Which  planets  would  be  suitable?  What  would  be  the 
advantages  and  disadvantages  of  tracking  a  planet  instead  of  a 
star? 


Answer :  Planets  can  be  used,  in  principle.  Most  suitable 

planets  are  Jupiter,  Saturn,  and  Venus  (when  its  angular  dis¬ 
tance  from  the  sun  is  great  enough).  The  main  advantage  is 
that  the  planet  is  a  very  bright  and  easily  identified  target. 
Disadvantages  include  geometric  limitations  (the  bright  planets 
all  lie  very  nearly  in  the  same  plane)  and  the  computational 
burden  associated  with  determining  the  planet's  position  as  a 
function  of  time. 


44 


Chapter  Five 


9.  Distinguish  between  active  and  passive  systems  for  tracking 
satellites  and  determining  their  orbits. 


Answer :  The  active  systems  transmit  energy  (radio,  radar,  laser) 

to  the  satellite,  resulting  in  reflection  or  retransmission  back 
to  the  ground.  Passive  systems  do  not  transmit,  but  receive 
energy  transmitted  by  the  satellite  or  reflected  from  it 
(as  in  the  case  of  sunlight).  When  the  terms  active  and  passive 
are  applied  to  the  satellite,  rather  than  the  tracking  system, 
the  distinction  is  between  satellites  that  have  onboard 
transmitters  and  those  that  do  not. 


10.  Discuss  the  similarities  and  differences  between  laser 
and  radar  systems. 

Answer :  In  principle,  they  are  the  same,  differing  only  in 

the  wavelength  of  the  energy  used.  A  laser  is  a  radar  operating 
in  the  visible  (or  infrared),  rather  than  at  so-called  radar 
frequencies.  There  are  practical  distinctions  in  terms  of  power 
supply,  instrumentation,  and  mounting  requirements.  Also,  lasers 
are  usually  restricted  to  clear-weather  operation. 

11.  Reference  is  made,  in  this  and  other  chapters,  to  radar 
devices  operating  at  C-band ,  L-band ,  S-band ,  etc.  To  what 
do  these  letters  refer? 

Answer :  An  old  designation  for  radar  frequencies,  going  back 

to  World  War  Two  security  requirements,  used  code  letters  like 
C-band,  L-band,  etc.  Although  this  system  has  been  officially 
obsolete  for  many  years,  it  is  still  used.  An  overview  of  t he 
sy  s  t  ein : 


BAND 

FREQUENCY  LIMITS 
(GHz  ) 

WAVELENGTH  LIMITS 
(cm) 

P 

0.3  to  1.0 

30.0  to  100.0 

L 

1.0  to  2.0 

15.0  to  30.0 

S 

2.0  to  4.0 

7.5  to  15.0 

C 

4.0  to  8.0 

3.75  to  7.5 

X 

8.0  to  12.5 

2.4  to  3.75 

K 

12.5  to  18.0 

1.67  to  2.4 

KU 

18.0  to  26.5 

1.1  to  1.67 

K 

a 

26.5  to  40.0 

0.75  to  1.1 

12.  The  orbital  inclinations  of  the  GEOS-3  and  SEASAT-1  satel¬ 
lites  are  given  in  the  text  as  115  deg  and  108  deg,  respectively. 
What  is  the  meaning  of  an  inclination  greater  than  90  deg? 

Answer :  Inclinations  greater  than  90  deg  indicate  a  retrograde 
orbit.  This  means  that  the  motion  of  the  satellite,  as  viewed 
by  a  hypothetical  observer  located  above  the  north  pole,  is 
clockwise ■  If  satellite  motion  is  counterclockwi se  (in  the 
same  sense  as  the  rotation  of  the  earth)  the  orbit  is  called 
direc t .  (The  term  prograde  is  sometimes  seen  instead  of  direct , 
but  this  is  an  incorrect  usage  and  should  be  avoided.) 

13.  Use  Eq.  (4.5-1)  to  estimate  the  velocity  of  an  ocean  current 
crossed  at  right  angles  by  the  subtrack  of  a  satellite  equipped 
with  a  radar  altimeter,  given  that  the  altimeter  detected  a 
change  in  ocean  height  of  1  id  over  a  300  km  current  width. 

Assume  that  the  crossing  occurred  at  latitude  30  deg  South. 

Answer :  Using  these  approximate  values  for  the  physical  constants 

2 

g  =  9.8  m/sec 

uj  =  7 . 3  x  10  5  rad/sec 
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and  the  given  data 
Ah  =  1  m 
4>  =  30  deg 
a  =  90  deg 
AL  =  3  x  105  m 

the  use  of  Eq .  (4.5-1)  gives 

V  =  0.45  m/sec 

Velocities  of  ocean  currents  are  often  expressed  in  knots. 
The  conversion  is 


1  knot  =  0.5144  m/sec 

Thus,  the  computed  velocity  may  also  be  given  as 

V  =  0.87  knot 

14.  Explain  the  meaning  of  the  word  ephemeris ■ 

Answer:  The  original  meaning  of  the  word  ephemeris  (plural 

ephemerides)  was  much  the  same  as  diary  or  journal ,  referring 
to  the  recording  of  information  on  a  daily  basis.  Astronomers 
used  the  word  for  tabulations,  on  a  daily  basis,  of  the  positions 
of  sun,  moon,  planets,  and  stars,  as  well  as  other  time-varying 
astronomical  data.  Current  usage  no  longer  requires  that  the 
tabulation  be  daily,  but  applies  to  the  representation  of  time- 
varying  information  in  general,  whether  by  listing  or  by  formula. 
In  the  special  sense  of  a  satellite  ephemeris ,  the  meaning  is 
a  representation  (by  formula,  tabulation,  or  otherwise)  of  the 
position  (often,  position  and  velocity)  components  as  functions 
of  time. 
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mm 


15.  Using  the  data  given  in  Table 
eccentricity  for  each  of  the  three 


A. 5-1,  calculate  the  orbital 
satellite  types  listed. 


Answer:  The  given  data  are  not  su 

since  the  latitude  of  perigee  and 
result,  the  proper  value  of  R,  the 
for  perigee  and  apogee  altitudes: 


fficient  for  a  precise  answer, 
apogee  is  not  given.  As  a 
earth’s  radius,  in  the  formula 


h  =  a(l+e)  *  R 
a 

h  =a(l-e)-R 
P 


is  not  known.  In  these 

e  is  the  eccentricity. 


formulas,  a  is  the  semi-major  axis  and 
Using  an  average  value  for  R, 


R 


6.37  *  106 


m 


(technically,  this  is 
volume  as  the  earth)  , 
are : 


the  radius  of  a  sphere  having  the  same 
and  the  data  from  Table  A. 5-1.  the  results 


— — — - T 

SATELLITE  TYPE 

eccentricity 

Navy  Navigation 

0.00 

High  Eccentricity 

0.75 

Low  Perigee 

0.37 

The  computation  is  done  most  easily  by  using  the 
in  the  form: 


formulas 


C  VhP 
a  -  R  +  2 


h  ->1 
a  p 

e  . 
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Chapter  Six 


16.  If  a  satellite  is  transmitting  at  a  frequency  held  precisely 
at  150  MHz,  and  is  approaching  an  observer  at  100  m/sec,  find 

the  frequency  of  the  signal  received  by  the  observer. 

Answer :  The  Doppler  shift  increases  the  received  frequency  by 

50  Hz  (Eq.  A. 6-2);  the  signal  as  received  by  the  observer  has 
a  frequency  of  150,000,050  Hz. 

17.  How  accurately  must  a  ground  observer  measure  the  frequency 
of  the  signal  received  from  a  satellite  in  order  to  determine 
the  radial  velocity  to  within  one  mm/sec?  Assume  that  the 
satellite  transmitter  frequency  is  A00  MHz. 

Answer :  From  Eq .  (A. 6-2)  it  is  seen  that  the  relation  between 

a  small  error  in  the  received  frequency,  6fr,  and  a  small  error 
in  the  radial  velocity,  6s,  is  given  by 


With  the  given  values 

ft 

6s 


A  x  108  Hz 
10  8  m/sec 


and  using 


c  =  3  x  108  m/sec 

for  the  velocity  of  light,  the  result  is 

6  f  =  0.0013  Hz 
r 
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This  level  of  accuracy  would  be  very  difficult  to  achieve  in 
practice . 

18.  During  a  tracking  interval  of  10  sec,  the  frequency  of 
the  signal  received  from  a  satellite  increases  linearly  from 

fr(0)  =  f re f  =  150  Mhz 
to 

£r{10)  =  f ref  +  50  Hz 

Find  the  total  Doppler  count  over  this  interval. 

Answer :  Using  Eq .  <4.6-4),  the  integrand,  f^-f  can  be 

expressed  as 


where  t,  the  time  measured  from  the  beginning  of  the  tracking 
period,  varies  from  0  to  10  sec.  Thus,  Eq .  (4.6-4)  becomes 

f10 

N12  =J0  dt 

and  N-p,  the  Doppler  count,  is  250  cycles. 

19.  Assume  that  a  ground  receiver  measures  a  Doppler  count  of 
250  cycles  during  a  tracking  interval  of  10  sec.  The  receiver's 
stabilized  reference  frequency  is  exactly  150  MHz,  while  the 
satellite's  transmission  frequency  is  known  to  be  149.998  MHz. 
Compute  the  change  in  range  (from  satellite  to  receiver)  during 
the  tracking  interval. 

Answer :  From  Eq .  (4.6-9),  the  change  in  range  is  computed  to 

be 

6s  =  -0.1  mm 
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For  all  practical  purposes,  the  change  in  range  is  zero.  This 
result  illustrates  the  fact  that  the  measured  Doppler  count 
has  two  sources: 

•  The  Doppler  shift  caused  by  relative  motion 

•  The  difference  between  transmitter  and  reference 
frequencies  (6f  =  f  -  fref)- 

In  this  example,  the  observed  Doppler  count  arises  almost  entirely 
from  the  difference  in  frequencies,  6f. 

20.  Discuss  the  accuracy  requirements  for  measured  Doppler 
count  in  terms  of  desired  accuracy  in  the  determination  of 
change  in  range. 

Answer:  Equation  (4.6-9)  shows  that  small  errors  in  the  change 

in  range ,  6s,  are  related  to  errors  in  the  Doppler  count,  N^. 
by  the  factor  c/f ref  For  the  NNSS  frequencies,  this  factor 
is  roughly  unity;  thus,  errors  in  N-^  (measured  in  cycles ) 
translate  directly  into  errors  in  6s  (measured  in  meters ) .  For 
example,  an  accuracy  level  of  0.01  cycle  corresponds  to  1  cm 
of  range  difference. 

21.  Why  does  the  NNSS  use  the  dual - frequency  (150  and  400 
MHz)  mode  of  operation? 

Answer :  The  use  of  two  frequencies  permits  compensation  for 

shifts  in  frequency  caused  by  ionospheric  effects.  Without 
this  compensation,  the  accuracy  of  the  system  would  be  degraded 
to  an  unacceptable  level. 


Chapter  Seven 


22.  Why  is  GPS  sometimes  referred  to  as  a  four -dimensional 
positioning  system? 

Answer :  The  terminology  reflects  the  fact  that  GPS  provides 

the  user  with  three-dimensional  position  information,  as  well 
as  precise  time . 

23.  Explain  the  meaning  of  the  term  pseudo-range ,  as  used  in 
the  context  of  GPS. 

Answer :  GPS  users  derive  range  information  from  measured  time- 

of-arrival  (TOA)  information.  Because  of  uncertainty  in  the 
user's  knowledge  of  GPS  system  time,  there  is  a  bias  error  in 
the  derived  range.  The  range  information  which  contains  this 
bias  error  is  called  pseudo-range . 

24.  Would  the  predicted  accuracy  of  the  fully  operational  GPS 
be  adequate  for  use  by  military  aircraft  as  a  precision  landing 
system  under  conditions  of  zero  ceiling  and  visibility? 

Answer :  No.  This  is  not  a  projected  application  of  GPS. 


Chapter  Eight 

25.  Why  are  remote-sensing  satellites  like  LANDSAT  placed  in 
nearly  polar  orbits? 

Answer :  To  permit  nearly  global  geographic  coverage.  For 

another  important  reason,  review  Exercises  4,  5,  and  6. 

26.  Which  agencies  of  the  United  States  government  are  primarily 
involved  in  the  processing  and  dissemination  of  LANDSAT  data? 

Answer :  The  U.S.  Geological  Survey  (under  the  Department  of 

the  Interior!  and  NASA  --  primarily  the  Goddard  Space  Flight 
Center  --share  responsibility  for  the  processing  and  dissemina¬ 
tion  of  LANDSAT  data.  Design,  launch,  and  control  are  a  NASA 
responsibility . 

27.  Identify  the  three  map  projections  used  for  the  presenta¬ 
tion  of  LANDSAT  data  in  map-compatible  form. 

Answer :  The  polar  stereographic  projection  is  preferred  for 

latitudes  higher  than  65  deg.  The  projections  used  at  lower 
latitudes  are  the  Universal  Transverse  Mercator  (UTM)  and  the 
Space  Oblique  Mercator  (SOM).  The  first  two  of  these  are  dis¬ 
cussed  in  Section  1.2.6  of  Unit  One.  A  detailed  discussion  of 
the  SOM  is  beyond  the  scope  and  level  of  these  Lecture  Notes; 
interested  students  may  be  referred  to  the  original  sources; 

•  "Space  Oblique  Mercator"  by  A.P.  Colvocoresses , 
published  in  Photogrammetric  Engineering, 

Volume  40  (1974),  pages  921  to  926. 

•  "The  Space  Oblique  Mercator  Projection"  by 
J.P.  Snyder,  published  in  Photogrammet ri c 
Engineering  and  Remote  Sensing,  Volume  44 

( ] 9 78 )  ,  pages  585- 596 . 


53 


The  following  simplified  explanation  may  be  helpful.  The  ordinary 
Mercator  Projection  involves  a  projection  from  the  globe  onto 
a  cylinder  whose  axis  coincides  with  the  earth's  spin  axis.  Globe 
and  cylinder  are  tangent  along  the  equator.  In  a  Transverse  Mercator 
Projection,  the  cylinder  is  tilted  90  deg,  so  that  its  axis  is 
perpendicular  to  the  earth's  spin  axis.  Globe  and  cylinder  are 
tangent  along  a  meridian.  For  the  Oblique  Mercator  Projection, 
the  cylinder  is  tilted  at  an  intermediate  angle  and  is  tangent 
along  some  great  circle.  If  this  great  circle  is  aligned  with 
the  subtrack  of  a  satellite,  the  resulting  projection  is  a 
first  approximation  to  the  Space  Oblique  Mercator.  Further 
stages  of  mathematical  refinement  are  required  because  the  actual 
satellite  subtrack  is  not  exactly  a  great  circle.  The  deviations 
arise  from 

•  Rotation  of  the  earth 

•  Precession  of  the  satellite's  orbital  plane 

•  Other  perturbations  in  the  motion  of  the 
satellite. 
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